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Abstract 



The control of thermal decoherence via dynamical decoupling and via the 
quantum Zeno effect (Zeno control) is investigated for a model of trapped ion, 
where the dynamics of two low lying hyperfine states undergoes decoherence 
due to the thermal interaction with an excited state. Dynamical decoupling 
is a procedure that consists in periodically driving the excited state, while the 
Zeno control consists in frequently measuring it. When the control frequency 
is high enough, decoherence is shown to be suppressed. Otherwise, both 
controls may accelerate decoherence. 



I. INTRODUCTION 



The theory of quantum information and computation has provided various promising 
ideas such as substantially faster algorithms than their classical counterparts and very se- 
cure cryptography [1]. Examples are Shor's factorizing algorithm [2] and Grover's search 
algorithm [3], where several computational states are simultaneously described by a sin- 
gle wave function and parallel information processing is carried out by unitary operations. 
Moreover, some of the basic steps have already been experimentally realized: Basic opera- 
tions for quantum computation were realized with trapped ions [4,5] and with the nuclear 
spins of organic molecules [6]. Shor's algorithm for factorizing Z — 15 was investigated with 
the nuclear spins of organic molecules [7]. 

The essential ingredient for the efficiency of quantum algorithms and cryptography is the 
principle of superposition of states. As pointed out e.g., by Unruh [8], the loss of purity (i.e., 
decoherence) of states would deteriorate the performance, particularly in the case of large 
scale computations or of long-distance communications. Thus, the information carried by a 
quantum system has to be protected from decoherence. So far, several schemes have been 
proposed, such as the use of quantum error-correcting codes [9], the use of decoherence free 
subspaces and/or noiseless subsystems [10] and the quantum dynamical decoupling [11-15]. 

The quantum dynamical decoupling was proposed by Viola and Lloyd [11], where the 
system is periodically driven with period Tc in an appropriate manner so that the target 
subsystem is decoupled from the environment. It was shown [11,12] that a complete decou- 
pling is achieved in the Tc ^ limit, or the limit of infinitely fast control. The procedure 
is simpler than the other methods because one only has to periodically drive the system. 
However, as it is not possible to achieve the Tc ^ limit, its performance for nonvanish- 
ing Tc should be investigated. Such studies were carried out for a two-level system in an 
environment via a system-energy- preserving interaction [11] and for a harmonic oscillator 
coupled with an environment [14]. Here, we will provide one more example, namely a model 
of a trapped ion used in Ref . [4] . This model explicitly involves a unitary operation for the 
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quantum-state manipulation, which was not included in the previous models. 

The key ingredient of dynamical decoupling is the continuous disturbance of the system, 
which suppresses the system-environment interaction. As already pointed out by Viola and 
Lloyd [11], the situation is similar to the so-called quantum Zeno effect, where frequent mea- 
surements of a system suppress quantum transitions [16-18] (for recent reviews, see [19]). 
This phenomenon is more general than originally thought: a nontrivial time evolution may 
occur in the case of frequent measurements under an appropriate setting. Namely, when 
the measurement process is described by a multidimensional projection operator, frequent 
measurements restrict the evolution within each subspace specified by the projection op- 
erator and a superselection rule dynamically arises [20]. Therefore, if one can design the 
measurement process so that different superselection sectors (defined by the given measure- 
ments) are coupled by the interaction between a target system and the environment, the 
system-environment interaction can be suppressed by frequent measurements. We refer to 
such a decoherence control as a quantum Zeno control. Since, in case of the quantum Zeno 
experiment by Itano et al. [17,18], the measurement process was realized as a dynamical 
process, namely the optical pulse irradiation, it is interesting to compare the two procedures 
(the quantum Zeno control and the quantum dynamical decoupling) for a model of trapped 
ion. This is one of the objectives of this article. 

This article will be organized as follows. In Sec. II, the quantum dynamical decoupling 
and the quantum Zeno control are briefly reviewed. In Sec. Ill, we introduce a model of the 
trapped ion, which takes into account the unitary Rabi oscillation and thermal decoherence. 
The dynamical decoupling and Zeno controls of this model are discussed, respectively, in 
Sees. IV and V. After discussing the cases of inflnitely fast controls, the effects of the 
finiteness of the control period are investigated. It is shown that both controls may accelerate 
decoherence if they are not sufficiently fast. This implies the necessity of a careful design 
of the control and a careful study of the timescales involved. The last section is devoted to 
the summary and discussion. 
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II. QUANTUM DYNAMICAL DECOUPLING AND QUANTUM ZENO 

CONTROL 



A. System 



The total system consists of a target system and a reservoir and its Hilbert space TYtot 
is the tensor product of the system Hilbert space, Tisi the reservoir Hilbert space, TIb'- 
Ti.^.^^ — Jig (g) TYb- The total Hamiltonian i^tot is the sum of the system part Hg ® 1b, the 
reservoir part I5 ® Hb and their interaction Hsb, which is responsible for decoherence: 



The operators I5 and Is are the identity operators, respectively, in the Hilbert spaces li.s 
and I-Lb-i and the operators Hg and Hb act, respectively, on Tis and TIb- Here, in order to 
discuss controls in an interaction picture, a time-dependent interaction is considered. 

Since, in general, the reservoir state is mixed, it is convenient to describe the time 
evolution in terms of density matrices. In the case of a quantum state manipulation, the 
initial state p(0) is set to be a tensor product of a system initial state (t(0) and a reservoir 
(usually equilibrium) state pB'- p(0) = (t(0) ®Pb- The system state (T{t) at time t is given by 
the partial trace of the state pit) of the whole system with respect to the reservoir degrees 
of freedom: a{t) = trBp(t). When a{t) is not unitarily equivalent to a{0) for a given class of 
initial states, decoherence is said to appear. The purpose of the control is to suppress such 
decoherence. For the decoherence control, it is sufficient to consider only those initial states 
which are relevant to the quantum state manipulation in question, but not all states. 



Here we slightly generalize the arguments of Ref. [12] (see also [14]). This control is 
carried out via a time dependent system Hamiltonian Hc{t): 



Htot ^ Hs®1b + 1s®Hb + HsBit) . 




B. Quantum Dynamical Decoupling 



(2) 
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where Hc(t) is designed so that Uc{t) = Texp |— i Jq Hc(s)ds^ satisfies 

(A) Uc{t) is periodic with period T^, Uc{t + Tc) = Uc{t). 

(B) J^^ dtiUlit) ® lB)HsB{t + s){Uc{t) ® 1b) = 0{T^+') , (0 < e < 1, : small, Vs). 
Going to the interaction picture where Hc{t) is unperturbed, the density matrix at time 

T = NTc with an initial state p(0) is given by p(T) = Utot{NTc)p{0)Uli{NTc) where 



f/tot(iVT,) = r exp -W H,,,{s)ds = n 



AT 



m=l 



Texp <j / Hiot{s)ds 

(m-l)Tc 



(3) 



and Htot{t) = {U^{t) lB)Htot{Uc{t) (8) Is). A standard Magnus expansion of the time 
ordered exponential [23] leads to 



Texp 



[ J{m-l)T, 



(0) , fr{l) 



where = ^ I(^i)Tc Htot{s)ds and the rest terms H^^ are of order of {j 
By assumption (B), one has H^'^ = H^^^ + 0{T^) where 



(4) 



1,2,---)- 



(5) 



Hs = ^ Jq" dtUl{t)HsUc{t) and they are independent of Tc because Uc{t) is Tc-periodic. 
Therefore, in the limit Tc — > while keeping T = A^Tc constant, one obtains 



C/tot(T) 



N 



(6) 



In short, as a result of the infinitely fast control, the system-reservoir coupling is eliminated 
and, thus, decoherence is suppressed. Note that if one designs He so that 



Ucit) = gj 



'-Tc<t < — Tc ;? = 1, 



M 



(7) 



where {gj} is a set of unitary operators acting on Tis, Hs becomes 



^ M 

Hs = j^Y.9]Hsgj 



(8) 



The relation between the dynamical decoupling with the prescription (7) and a symmetry 
group was discussed in Ref. [13]. 



C. Quantum Zeno control 



Now we turn to the Zeno control by adapting the argument of Ref. [20]. This control 
is performed by frequent measurements of the system. The most general measurement is 
described by a projection operator acting on the density matrix: 



where {Pn} is a set of orthogonal Hermitian projection operators acting on Tis- In the 
following, we restrict ourselves to the case where the measuring apparatus does not "se- 
lect" different outcomes (nonselective measurement) [21] and the projection operators are 
complete; J2n Pn = Is- As in the dynamical decoupling, the measurement is designed so 
that 



The Zeno control consists in performing repeated nonselective measurements at times 
t — nTc {n = 0, 1, 2, • • •). Between successive measurements, the system evolves via i^tot- In 
terms of the Liouville operator Ctot defined by AotP = [-f^tot, p] = -f^totP — pH tot, the density 
matrix p{NTc) after + 1 measurements with an initial state p(0) is given by 



Pp = Y,{Pn® U) P {Pn ® U) . 



(9) 



n 



(C) PHsB{t) = En (^n ® 1b) HsB{t) (P„ ® 1b) = 0. 




m=l 



(10) 



Assumption (C) yields 




nP[H's 1b + I5 Hb, p] = nt'^p , 



and, thus, in the hmit Tc — > while keeping T — NTc constant, we get 



p{NTc) ~ p n {(^-'p 



[ 




N 



P(0) 



m=l I 




(11) 



where the controlled system Hamiltonian H'c^ 



'g is given by 
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H's = J2PnHsPn. (12) 

n 

Hence, as a result of infinitely frequent measurements, the system-reservoir coupling is elim- 
inated and, thus, decoherence is suppressed. Note the similarity between the controlled 
system Hamiltonians for a particular dynamical decoupling, (8), and for the Zeno control, 
(12). This is not a mere coincidence. Indeed, one can show that, by enlarging the Hilbert 
space so that the original measurement process is expressed by a dynamical process in the 
larger space, the two controls are equivalent. This will be discussed in detail elsewhere. 
However, throughout this article, the dynamical decoupling refers to a situation where the 
evolution is coherent (unitary) and the Zeno control to a situation where the evolution 
involves incoherent processes such as measurements. 

III. A TWO-LEVEL SYSTEM WITH THERMAL DECOHERENCE 

A. Model 

We consider the model of a trapped Be ion used in Ref. [4] (see also [17]). Here we assume 
that the ion is at rest and consider only the dynamics of the hyper fine states. Then, the 
main mechanism of decoherence can be attributed to the emission and absorption of thermal 
photons associated with transitions to nearby excited states. For the sake of simplicity, one 
of the hyperfine states is assumed to couple electromagnetically with a nearby excited state 
(see Fig. 1), the polarization of the photon is neglected and only the rotating terms (i.e., 
the slowly varying terms in the interaction picture) are taken into account (rotating wave 
approximation) for the driven parts. 

Let |1), 1 2) and |3) be the lower hyperfine, upper hyperfine and excited states, respec- 
tively, and Ok the annihilation operator of a photon with wave vector k and energy Uk = c|k| 
(c: speed of light). Then, the Hamiltonian is given by 
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^\/\y\/v 

thermal \ y > |2> 

P^°'°" Rabi Oscillation A 

FIG. 1. A schematic picture of the system. |1), |2) and |3) are the lower hyperfine, upper 
hyperfine and excited states, respectively. 

i7tot(i) = iis{t) (8) 1b + I5 + \iisB , (13) 

Es{t) = u;2\2){2\ + cc;3|3)(3| + A2{Ae*-^*|i)(2| + (h.c.)} , (14) 

Hb ^ j d^k ujk ajk^k , (15) 
HsB = Jd'm (|1)(3| + |3)(1|)® (aj, + ak) . (16) 

The bare energies 0J2 and 103 of states |2) and |3) are measured from that of the lower 
hyperfine state |1) {uui = 0). The third term of (14) represents the rf-control of the Rabi 
oscillation and the amphtude A^A corresponds to the Rabi frequency. Because of the Lamb 
shift, the energy difference between the two hyperfine states is different from 002 and the 
frequency of the irradiated field should be so tuned that it resonates with lv'2- The function 
Vk is assumed to behave like 



with a cut-off frequency lUc and a dimensionless strength vq. The dimensionless coupling 
constant A measures the relative order of magnitude of each term and is of order ^^g/us, 
where 7e is the inverse lifetime of the excited state. For the system in [4] , typical order of 
magnitudes of the frequencies are 7e/co'3 ~ 10~^ and (Rabi frequency) /a;3 ~ 10~^°, which 
imply that A ~ 10~^ and the Rabi frequency is of order A^. Also it is reported that the 
decoherence time is longer than the Rabi period by two orders of magnitude [4] . 
It is convenient to move to a rotating frame with the aid of the unitary operator 

UR{t) = exp {i (^X^d\l){l\+cu2\2){2\+uJs\3){3\)t}®exp{iuj'^t j d^k alaA , 



where — UJ2 — 0J2 and uj'.^ — uj^ — \^5. Then the transformed Hamiltonian is 



H^^-S\l){l\ + {A\l){2\ + {h.c.)}, (19) 

I (fk {iUk - 4) 4«k , (20) 

H^Bit) = / d^m |1)(3| (at + e-^"^'^W) + (h.c.) (21) 

This is our starting point. 



B. Decoherence 

We consider the time evolution starting from an initial state given by the tensor product 
of a system initial state and the reservoir equilibrium state with inverse temperature /3: 



= {Cb + XCsBit) + X'Cs} Pit) (22) 
p(0) = a(0) ® pB (23) 
Pb^^bM-PHb) (24) 

where Z is the normalization constant and the operators Cb, Csb and Cs are defined by 

CbP = [Is ® Hi p] , CsB{t)p = [HiB{t),p] , Csp = [Hi ® 1^, p] . (25) 



Since the time scale of the quantum state manipulation is of the same order of magnitude 
as the Rabi period (~ A~^) and is very long compared with l/cu^ (~ A°), the process 
is well described by the van Hove hmit approximation [24,25]. The starting point is the 
decomposition of the Liouville equation (22) with the aid of a projection operator 

Vp = {tiBp) ® pB (26) 

where tr^ stands for the partial trace over the reservoir degrees of freedom and pb is the 
equilibrium reservoir state (24). Then, in the limit A ^ while keeping r = X^t constant, 
the reduced density matrix o" = tr^p is found to satisfy [24-26] 
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^ = -^Aa , (27) 

where we have used VCbV — and 

-iKa = a\ - dt lim trs \Csb (-^) e-'^^*£5B - e'^«* (7(r) psj 

= -z[5'|3)(3| + A(|l)(2| + (M),<7] 

+7e|l)(3|a|3)(l| + 7d|3)(l|<7|l)(3|-{||l)(l| + ||3)(3|,a} , (28) 

with { , } the anti-commutator. In the above, the hmit A — > should be understood to drop 
terms which oscillate with frequencies ~ (cf. [25]). The parameter 5 is chosen as 

S^V.v.Tdu^ (29) 

and S', 7d and are given by 

S' = -p.v. r .1.-^ (30) 

J—oo CO — Cc'3 

7d = 27r«;d('^3)' 7e = 27rfi;e('^3) > (31) 

where 

^ ^ J k ^pu,k _ 1 27r - 1 ^ ^ 

are the thermal spectral density functions (form factors). The symbol p. v. in front of the 
integrals indicates Cauchy's principal value. 

In terms of the matrix elements (7^ = (i|(7|j), one has (T21 = ^12 and 

= -iA{(72l - (712} - 7<i(7ii + 7e(733 (33) 
= -lA{a22 - CTu} - yCri2 (34) 

da22 



dr 

9(733 



= iA{(72i - (712} (35) 
= 1d(7ll - 7eCr33 ■ (36) 



dr 

The purity of the target states is measured hj rj = afi + 0"|2 + 2|cri2p, as = 1 for pure 
superpositions of |1) and |2), and 77 < 1 for states involving the irrelevant state |3) or mixed 
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states. In Fig. 2, the evolution of the quantity r]{t) starting from a{0) — is shown 

for A — 1007rf, 7e = 10007^. As time goes on, the purity 77 of the target states is lost, or 
decoherence takes place. One clearly sees that the decoherence time scales as 1/7^. 



0.8 
0.6 
0.4 
0.2 


2 4 6 8 10 

t 

FIG. 2. Time evolution of the purity rj of the target states. The time unit on the horizontal 
axis is the decoherence time 7^^. 




IV. CONTROL OF THERMAL DECOHERENCE VIA DYNAMICAL 

DECOUPLING 

We consider a dynamical decoupling control of the thermal decoherence discussed in the 
previous section. Since decoherence arises from the transition between the states |1) and |3) 
associated with absorption and emission of photons, it is expected to be suppressed if |3) 
does not contribute to the |1)-|2) dynamics. So we consider a control via the Rabi oscillation 
between the state |3) and a higher excited state |4) (cf. Fig. 3), 

14> 

Rabi Oscillation O. 

|3> 




^\/V/\/V' 

thermal \ , , |2> 

photon \ / J Rabi Oscillation A 



FIG. 3. A schematic picture of the system under the quantum dynamical decouphng control. 
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which is described by 

H^{t) = Qe^('"4+«)*|3)(4| + (h.c.) , (37) 

where 004 is the energy of the state |4), and Q and ^ are real parameters. Since there exists 
one more state, a term c<j4|4)(4| ® 1b should be added to the Hamiltonian Hs- As a result, 
the following Liouvillian should be added 

{ACs + Cc{t)}p=[{uj4A){A\ + H,{t))®lB , p] , (38) 

and the evolution equation reads 

^ = {Cb + ACs + C,{t) + XCsBit) + X'Cs} p{t) . (39) 

A. Ideal dynamical decoupling 

The evolution operator Uc{t) generated by the control Hamiltonian Hc{t) is given by 

U^{t) = r exp (^-i dt'H^{t')^ 

= e-'('^^+«l^><^l* expH{-(a;4 + 0|4)(4| + fi(|3)(4| + (/i.c))}t] . (40) 

When restricting to a subspace spanned by |3) and |4), the second factor is a sum of two 
oscillating projection operators with frequencies 

a;± = ^ (-004 - e ± + 0^ + 4^^') (41) 

and the first factor is a sum of a time-independent and an oscillating projection operators, 
the frequency of the latter being 104 + ^. Thus, Uc{t) is a sum of four oscillating terms with 
frequencies ±a;+, ±u;_ and is Tp-periodic provided both uj+Tc and lv-Tc are integer multiples 
of 27r. 

Under this prescription, {U}.{t) (g) \B)HsB{Uc{t) ®1b) is a sum of terms proportional to 
g±ia)+t^ e^^'^-*, e±2i^te±^'^+*, e^'^^'^^e^'^-^ and its average vanishes in the ^ hmit: 
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lim 1 dt{Ul{t) ® lB)HsB{Uc{t) ® 1b) = . (42) 

Therefore, the general argument of Sec. II shows that the coupUng between the system and 
the reservoir is suppressed in the limit Tc — > (that is a;± — > oo) while keeping T = NTc 
constant. Moreover, the system obeys the Hamiltonian 

= ^ {^^1^)^^! + ^'^s) Uc{t) 

= ^^l.^^s (Q2 + ^2)2 + >^Hs . (43) 

Therefore, the target system spanned by |1) and |2) is free from decoherence and performs 
ideal Rabi oscillations. 

It is interesting to see the relation between the dynamical decouphng and a dynamical 
quantum Zeno effect due to a "continuous" measurement [20]. For the particular choice 
(37), it is possible to eliminate the explicit time dependence of the control Hamiltonian He, 
by going to another rotating frame with the aid of the unitary operator 

U'j,{t) = exp {t{uj, + 0|4)(4|t} ® 1b , (44) 
and, then, the transformed density matrix p — Uji{t)pU'^{t) obeys 

^ = {Cb + AC's + 4 + ^^SB + X'Cs} P (45) 
where the transformed control Liouvillian is 

{A4 + 4}p=[(-e|4)(4| + Q{|3)(4| + (/i.c.)})®ls , p\ . (46) 

In this picture, state |4) of energy — ^ is coupled by a constant couphng Q to state |3). The 
short-period limit Tc corresponds to the strong coupling limit fi, ^ — >• oo, because Cj±Tc 
must be integer multiples of 27r. But this is just the case of a dynamical quantum Zeno effect 
due to a continuous measurement [20], where i^meas = ~^|4)(4| + Jl{|3)(4| + (/i.e.)} plays 
the role of a measurement Hamiltonian. One can show that, in the hmit of strong coupling, 
a dynamical superselection rule arises, the Hilbert space is split into Zeno subspaces and 
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the system Hamiltonian is given again by (12), where the projections P„s, defining the Zeno 
subspaces, are the eigenprojections of i^meas [20]. This is a consequence of an interesting 
relation between strong-couphng regime and adiabatic evolution [22]. 
The eigenprojections of i^meas are given by P± — |±)(±|, where 

|±> = "^^^ , (47) 



with eigenvalues 



a;± = i (-e ± + ■ (48) 

Therefore from (12) one gets 



which is nothing but the Hamiltonian (43) under the transformation (44). We therefore see 
that, in this particular case, dynamical decoupling is completely equivalent to the dynamical 
Zeno effect. 



B. Nonideal dynamical decoupling 

Here we consider the case of nonvanishing Tc and solve the evolution equation (39). As 
pointed out in [11,14], the ideal dynamical decoupling is achieved when the control frequency 
27r/Tc is higher than the threshold frequency cUc in the system- reservoir interaction Hsb and, 
thus, we consider the case where 27r/Tc = 0(A°). Then, the slow process which is relevant 
to the quantum state manipulation is well described by the van Hove limit approximation 
[24-26]. 

We consider the evolution in the rotated frame (45). Note that, since the transformation 
Uji in (44) does not affect the evolution of the states \ j) {j = 1, 2) and the field variable, one 
has {i\trBp\j) — (i]trBp]j) — 1,2). By the standard procedure of the van Hove limit 
approximation [24,25], in the limit A — > while keeping r = XH constant, one obtains 
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^ = -iAB(T (50) 



where (t(t) = trsp and 



— iAR(7 = — i 



$:5,|s)(s|+A(|l)(2| + (/i.c)),a 



.s=± 



+ E(7l|l)(«W(l| + 7jk)(lHl>(«l)-|^|l)(l| + E|l«)(«l>'^| ■ (51) 

The states |±) are the normahzed hnear combinations of the states |3) and |4) given by (47) 
and the decay rates 7| and 7^ (s = ±) and 7^ are given by 

Id-ld+ll, 7je = 27r-^^«d/e(4 + ^±), (52) 

' CO^ — CO— 

where u!± are the frequencies (48) and Hd/eil^) are the thermal form factors (32), extended 
to the whole real axis due to the counter-rotating terms. (Incidentally, notice the exchange 
symmetry Kf,{uj) = Kd{—uj) of the extended form factors.) 

The prefactors in the second equation in (52) are nothing but the squares of the matrix 
elements between the undressed state |3) and the dressed states |±) (47): 

'"-I --±)\\ (53) 



a;+ — a;_ 

The explicit expressions of the Lamb shifts 5s {s — ±) of the excited states are omitted 
since the relevant sector of the evolution equation does not depend on them. Note that the 
parameter 5 is chosen so that the operator Ab does not contain a commutator with |1)(1|. 
In terms of the matrix elements (Tij = one has 

^ = -zA{a2i - a,,} - ^^an + E ^>ss (54) 

^'^ s=± 

^ = -tA{a22 - an} - (55) 

= iA{(T2i - (712} (56) 

^ = 7^^ii - 7eV.. . {s = ±) (57) 

The evolution of the purity rj = afi + cr|2 + 2|cri2p of the target states is shown in Fig. 4 for 
different values of control parameters, where the parameters in the control Hamiltonian are 
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set to ^ = 24Q/5 (which gives a;+ = Q/5 and uj- — — 5Q), the frequency-cutoff to ujc — 10a;3, 
and the other parameters are chosen so that one has A — 1007^ and 7e = 10007^ for the 
uncontrolled case. As in the previous section, the initial state is (t(0) = |1)(1|. Fig. 4 shows 
that the dynamical decoupling control may accelerate decoherence if the parameters are not 
appropriately tuned. 

1 

0.8 

0.6 
0.4 
0.2 


2 4 6 8 10 

t 

FIG. 4. Evolution of the purity r] of the system state. The time unit in the horizontal axis 
is the decoherence time for the uncontrolled case, (a) Control frequency |a;_| = 150 x w'^; 
(b) control frequency = 0.5 x uj'^. For comparison, the behavior of rj without control is also 
displayed (broken curve). 

In Fig. 5, the control-frequency dependence of the decoherence rate 7^ is shown. As can 
be seen in the figure, decoherence is first enhanced, for small values of = 5J1, and then 
suppressed for much larger values of |a;_|. Since the decoherence rate 7^" due to |-|-) is a 
monotonically decreasing function of this can be understood as follows: In the rotating 
frame where the |3)-|4) oscillation is ehminated, when Q = 0, state |3) is separated from 
the decay product state |1) by an energy 00'^. When f2 is turned on, state |3) splits into two 
dressed states |±) which are separated from |1) by the energies -|- a;+ = -|- |a;_|/25 and 
Lo'^ — \iO-\, respectively. The latter state is closer to state |1) than in the uncontrolled case 
and leads to a shorter decoherence time provided cug — |a;_| > 0. This is the deterioration 
observed in Fig. 4, case (b). On the other hand, if exceeds a threshold energy ujth = uj'^, 
the energy of state |— ) becomes lower than that of state In such a case, the counter- 
rotating term (which are now "rotating") does contribute to the decoherence rate. Notice 
that now, being uj'^, — \uj-\ < 0, 7" is smaller than 7^, as it should. Even after |a;_| has 
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exceeded the threshold cuth, ^^^^^ increases with \u!-\, since the state |1) is now unstable. 
And, finally, when the two dressed states are sufficiently far apart from level |1), the decay 
rates (and therefore decoherence) are suppressed because of the high energy cut-off of the 
form factor (17). Such values of |a;_| ~ 80a;3 is extremely higher than the threshold (Vth — u;^ 
(Fig. 5) and involve extremely short timescales [28]. 
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FIG. 5. Decoherence rate 7^ vs control frequency \uj-\/uj'^. 




V. QUANTUM ZENO CONTROL OF THERMAL DECOHERENCE 

For the same reason as in the dynamical decoupling, we disturb the evohition of |3) by 
repeated measurements (Fig. 6). The nonselective measurement of |3) causes the following 
change of the density matrix: 

p ^ Pp = 7r3p7r3 + (1 - 7r3)p(l - TTs) (58) 

where 713 is a projection operator acting on the whole Hilbert space tts = |3)(3| ® 1^ and 
1 = I5 <8) 1b- Then, the density matrix under the Zeno control is given by 

p{NT,) = n |pre"'CVc^'-W'^*pjp(0) , (59) 

m=l 

where Tc stands for the time interval between successive measurements. 




^\y\y\y\j' 

thermal \ , > |2> 

photon \ / ; Rabi Oscillation A 



17 



FIG. 6. A schematic picture of the system under the quantum Zeno control. 



A. Ideal Zeno control 

First we consider the case where Tc while keeping T = NTc constant. Then, as 
discussed in Sec.II and because PHsb — t^2,Hsbt^z + (1 ~ '^z)Hsb{'^ — tts) = 0, the state at 
time T is given by 

p{T) = P (e-'^'s'^ (8) e-^^-^p(0)e^^^^ ® e^^-^) , (60) 

where the controlled system Hamiltonian H'g is given by 

H's = {P^HiP^ + {Is - Pz)Hi{ls - P3)} 

= A2{-5|l)(l| + (A|l)(2| + (/..c.))} , (61) 

with P3 = |3)(3|. Hence, as a result of infinitely frequent measurements of state |3), the 
system-reservoir coupling is eliminated and, thus, decoherence is suppressed. 



B. Nonideal Zeno control 

As in the dynamical decoupling, we consider the case where Tc ~ l/tUg ~ A°, so that the 
time evolution is well described by the van Hove limit where A ^ while keeping r = X^NTc 
and Tc constant. We are here looking at the subtle effects on the decay rate arising from 
the presence of a short-time quadratic (Zeno) region. Therefore, it is important to notice 
that the standard method [24,25] is not applicable to the present situation and the limit is 
evaluated as follows: 

(1) Second order perturbation, up to A^, and PHsb — lead to 

/•(m+l)Tc rt , ^ ^ ^ 

-AM dt dse'^^'CsB{t)e-'^^^'-'^CsB(s)e-'^'''P\ . (62) 

(2) In terms of the operator K.^^, defined as a solution of the operator equation 
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^ r{m+l)Tc rt ^ ^ , „ ^ ^ r{m+l)Tc 

P dt dse^^^*£5B(t)e-'^^(*"")/:ss(s)e-'^«^P = P / dte*^«*/C^e-'^«* , (63) 

JmTc JmTc JmTc 

one has 

where /C(i) = /C„ for mi; < t < (m + l)Tc. 

(3) With the aid of (64) and pa = P^aP^ + (Is - ^3)^^(15 - P^), the final reduced state 
(t{t) is given by 

a(r) = lim tr^ i H |pre"'C'i)T. ^-W'^^pj ^(q) ® 

r=A2ArTc:finite \m=l ^ 'J 

= p hm trBp*(A^7;) . (65) 

r=A2iVTc:flnite 

where p\NT,) = pre-'C'(^^+^'^«-'^''=(*»''V(0) ® ps- 

(4) As p*(t) is a solution of 

- -i{Cb + A^^s - ^A2/C(^))p*(^) , p*(0) = a(0) ® pb , (66) 
the standard van Hove limit arguments [24,25] show that a(r) = pcr*(r) and a* satisfies 
. _i trs {(£5 - ^/C) a*(r) ® p^} , a*(0) = a(0) , (67) 



dr 

where the time-dependence of /C is lost as a result of the partial trace. 
As in the previous sections, the parameter 6 is chosen so that the |l)(l|-term does not 
appear in the evolution operator of a*. In terms of the matrix elements a*j = {i\a*\j), (67) 



reads 

da 



dr 



ii = -iA{a;, - al,} - ^lal, + 7f a*3 (68) 



''^^^ - -iA{al, - al,} - ^al, (69) 



dr ' 2 

tA{a;, - al,} (70) 



dr 

'-"-^33 rs,Z^* ('71 ^ 

-Q^ = Id (^11 - 7e C^33 ' ("IJ 

where the decoherence rate 7^ and the inverse lifetime 7^ of |3) are given by 
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7e 



7f = Tc / d(jj Kd{uj) sinc^ 



/■oo 



sine 



2 




(73) 



(72) 



where Kd/ei!^) again the (extended) thermal form factors (32) and sine(a;) = (sina;)/a;. 
The deeay rate 7^ in (72) should be compared to 7^ in (52). They express the (inverse) 
quantum Zeno effect, given by pulsed or continuous measurement, respectively [19]. 

Since the projection operator p does not affect the |1)-|2) sector, one has cr*^(r) = 
{i\a{T)\j) for a class of initial states where only the matrix elements {i\a{Qi)\j) = 1,2) 
are nonvanishing. Hence, rj — + + 2|cri2p measures the purity of the target states. 
Its evolution is shown in Fig. 7 for different values of 27r/Tc, where tOc — l^to'^ and the other 
parameters are chosen so that one has A = 1007^ and 7e = 10007^ for the uncontrolled case. 
As in the previous sections, the initial state is (t(0) = Fig- 7 shows that the Zeno 

control may accelerate decoherence if the parameters are not appropriately chosen. This 
can be seen more clearly in the control-frequency dependence of the decoherence rate 7^ , 
which is shown in Fig. 8. When the control frequency 27r/T'c belongs to a certain range, 
decoherence is enhanced. 

The enhancement of decoherence is qualitatively similar to the case of the dynamical 
decoupling. However, the high frequency behavior of the decoherence rate and its peak 
values are quite different. The high-frequency decoherence rates 7^ and 7^ , respectively, 
for the dynamical decoupling and Zeno control, are approximated by 



Therefore, 7^ decays exponentially for large |a;_| because of the exponential cut-off of the 
form factor and may take a maximum of order uj^^^^c/ ~ 140. On the other 
hand, 7f decays polynomially for large 2ti/Tc and 7^ could be much larger than 7^ because 
^^Jljoj'^ ~ 10^ is very large. 



la- 




e 





(74) 
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FIG. 7. Evolution of the purity r] of the target states. The time unit on the horizontal axis is 
the decoherence time for the uncontrolled case, (a) Control frequency 2-k/Tc = 5 x 10^ x uj'^; 
(b) control frequency 271 /Tc = 0.5 x io'^. For comparison, the behavior of rj without control is also 
displayed by a broken curve. 
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FIG. 8. Decoherence rate 7^ vs control frequency 21^ / {TcOo'^^) ■ 



VI. CONCLUSIONS 



In this article, we have studied the dynamical decoupling and Zeno controls for a model 
of trapped ions, where decoherence appears in the dynamics of the hyperfine states due 
to emission and absorption of thermal photons associated with the transition between the 
lower hyperfine and an excited state. By very rapidly driving or very frequently measuring 
the excited state, decoherence is shown to be suppressed. However, if the frequency of the 
controls are not high enough, the controls may accelerate the decoherence process and may 
deteriorate the performance of the quantum state manipulation. 
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The acceleration of decoherence is analogous to the inverse Zeno effect, namely the 
acceleration of the decay of an unstable state due to frequent measurements [27]. In the 
original discussion of the Zeno effect [16-19], very frequently repeated measurements of an 
unstable state is shown to slow down its decay. But, if the duration between two successive 
measurements is not short enough, the frequent measurements may accelerate the decay. 
This is the inverse Zeno effect. Obviously, this situation precisely corresponds to the increase 
of decoherence observed in this article. Moreover, since a very intense field is used for the 
dynamical decoupling control, the decrease of the decoherence time is also a consequence of 
the decrease of the lifetime of the unstable states due to the intense field [28] . 

There is room for improvement and further analysis: a number of neglected effects can 
be considered, such as the role of counter-rotating terms and Fano states, the influence of 
the other atomic states, the primary importance of the relevant timescales, and so on. These 
aspects will be discussed elsewhere. 
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